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PREFACE 


The SRI Automatic Tracking System (SATS) Is an objective cloud- 
tracking technique for automatically computing cloud motion vectors from 
digital picture data recorded by geosynchronous weather satellites. It 
Includes steps that select potential tracers from each frame, Identify 
the same tracer In successive f-ames, and convert the displacements into 
motion vectors. Since the computations are made without Intervention by 
the operator, the system has good potential for coping with large volumes 
of data for research and operational uses. 

In this study, we have significantly improved the automatic system. 

One major improvement provides better targets for tracking, particularly 
at cloud edges or In other areas of strong gradients In the data. A 
second change minimizes the number of spurious vectors computed near area 
boundaries. A third change allows the MOTION routine to calculate cloud 
displacements In as many as four separate layers simultaneously. The 
improvements were tested using data for complicated multilevel cloud 
motions. The cloud motions took place in the vicinity of the eye of 
hurricane Eloise as it passed across the Gulf of Mexico on September 22, 
1975. On this data, NASA Goddard personnel obtained rapid-scan observa- 
tions separated by 10 minutes or less. These images show that just to 
the northwest of the eye, there was high-speed outflow at the cirrus 
level. The cirrus shield had a sharp north-south edge. To the west of 
the cirrus edge, low clouds were moving In the opposite direction. 

Using infrared data for target selection, the automatic system was able 
to follow these motions. To the south of the eye, clouds were mainly 
cumuloform, and the automatic tracking following the motions using either 
visible or Infrared data in the target selection (grouping) process. 

The method was successful also in tracking motions in other nearby sectors 
of the hurricane. 

We have evaluated the automatic motions qualitatively by comparing 
them to the flow shown in a rapid- scan movie of Eloise made available to 
us by NASA-Goddard personnel. We have also compared the computer- 
determined motions qualitatively to those measured by NASA analysts and 
there is good overall agreement. For data with 4-km and 8-km resolution, 
the automatic system obtained a coverage of motions that was very similar 
to that obtained by the human analysts using data of comparable resolution. 

We have introduced a cross-correlation computation as an option in 
SATS and have compared cloud displacements computed by the MOTION program 
and by the cross-correiation technique. The cross-correlation computa- 
tion uses the visibltj or infrared targets selected automatically, and 
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finds their point of best fit by using a search routine. The accuracy 
of Che results is similar to that of the MOTION routine; however, cross 
correlation has the advantage that it gives a denser coverage of cloud 
motion vectors. 

We believe that the cloud motion vectors computed using the SRI 
automatic system are now competitive in accuracy and coverage with 
motion fields dot'o»'mined by human analysts working within reasonable 
time limits. We recommend that the automatic system be included as an 
option on AOIPS so that larger scale testing and evaluation can be 
carried out. 


CONTENTS 


PREFACE . ii 

LIST OF ILLUSTRATIONS ....... v 

I INTRODUCTION 1 

II CHANGES MADE TO THE AUTOMATIC METHOD 2 

III INTRODUCTION OF A CROSS-CORRELATION COMPUTATION .... 11 

IV FURTHER RESULTS OF COMPUTATIONS ..... 13 

V SUMMARY AND CONCLUSIONS 26 

REFERENCES 27 


iv 


ILLUSTRATIONS 


1 Block Diagram of the Automatic Cloud Tracking System .... 3 

2 Visible Image of the Northwest Sector of Hurricane 

Eloise at 1932 GMT, September 22, 1975 ........... 4 

3 Infrared Image for the Same Time and Area as Figure 2 ... 4 

4 Image of the Positive Infrared Deviations (Local Values 

Minus the Smoothed Background) 6 

5 Cloud Targets (Touching Groups) Based on Infrared Data 

After One Iteration of the DELETE Subroutine 6 

6 Cloud Targets After Two Iterations 7 

7 Cloud Targets After Three Iteraticns and With Groups 

That Abut Edges Deleted ..... ........ 7 

8 Cloud Motion Vectors Computed From the Groups of Figure 7 
and the Next Rapid-Scan Image Approximately 10 Minutes 

Later, Showing Cirrus Outflow and Opposite Low-Level Flow . 8 

9 Motion Vectors Approximately One-Half Hour After Figure 8 . 8 

10 Motion Vectors Computed Using Infrared Images Separated 

by 30 Minutes ............ ..... 9 

11 The Series of Cloud Motion Vectors Computed Using Six 
Successive Rapid-Scan Images Except for the Last Pair, 

Which Has the Normal 30-Minute Interval 9 

12 Cloud Motions Computed Using Visual Grouping of 8-km 

Resolution data for the Same Picture Pair as in Figure 8 . . 10 

13 Cloud Morions Computed by Cross Correlation for the Same 

Groups of Visible Data Used in Figure 12 ......... . 12 

14 Cloud Motions Computed for the Left Central Portion 
of Figure 3 (Including the Edge of the Cirrus Sheet) 

Using 4-km Infrared Data 14 

15 Cloud Motions Computed by Cross Correlation for the 

Same Infrared Groups Used in Figure 14 ... 14 

16 Cloud Motions for the Second Picture Pair Computed 

Using the MOTION Routine 15 

17 Cloud Motions Computed by Cross Correlation for the 

Same Infrared Groups Used in Figure 16 . 15 

18 Cloud Motions for the Third Picture Pair Computed 

Using the MOTION Routine .................. 16 


V 


19 Cloud Motions Computed by CrosH Correlation for the 

Same Infrared Groups Uaed in Figure 18 . 16 

20 Cloud Motions for the Series of Six Images of 4~km 

Infrared Data for the Left Central Portion of Figure 3, 

Computed Using the MOTION Routine , , 17 

21 Cloud Motions Computed by Cross Correlation for the Same 

Groups Used in Figure 20 ................. . 17 

22 Visible Image Using 4-km Data for an Area 280 by 480 km 

to the Soutliwest of the Eye of Kloise ........... 18 

23 Image of 4-km Infrared Data Co>*res,4)nding to Figure 22 . . . 18 

24 Cloud Motions Computed by the MOTION Routine Using 

Visible Targets for the First Picture Pair for the 

Area of Figure 22 19 

25 Cloud Motions Computed by Cross Correlation for the 

Same Groups Used in Figure 24 19 

26 Cloud Motions Computed Using Visible Targets for the 

Second Picture Pair .................... 20 

27 Cloud Motions Computed by Cross Correlation for the 

Same Groups Used in Figure 26 20 

28 Cloud Motions Computed Using Visible Targets for the 

Third Picture Pair 22 

29 Cloud Motions Computed by Cross Correlation for the 

Same Groups Used in Figure 28 ............... 22 

30 Sequence of Motions of Visible Targets for Five 

Picture Pairs, Including a 30-Minute Interval ....... 23 

31 Sequence of Motions for the Visible Targets of 

Figure 30, Computed by Cross Correlation ..... 23 

32 Sequence of Motions for Five Picture Pairs 

Computed by the MOTION Routine Using Infrared Data 24 

TABLES 

1 Comparison of Cloud Motions Computed Using the MOTION 

and CORRELATION Routines 23 


vi 


1 INTRODUCTION 


With the recent deployment a f four or more geosynchronous weather 
satellites around the globe, a very large amount of data concerning the 
earth's cloud cover is acquired each day. The data are n»re than ade- 
quate for achieving the long-term goal of determining cloud motions in 
real time over much of the globe, but the present processing methods are 
far too slow. This report describes our most recent work toward improv- 
ing the automatic cloud tracking methodology that has been under develop- 
ment for several years (Endlich et al., 1971; Wolf et al., 1977). We 
believe that an automatic method of this sort combined with recent 
advances in computers and display systems can significantly enhance the 
usefulness of geosynchronous satellite observations for research, and 
probably also for operations with consequent economic benefits. 

We have made significant improvements to the automatic system. The 
improvements are demonstrated in tests with a small sample of data for 
complex cloud motions in the vicinity of the eye of hurricane Elolse. 
Special rapid-scan data were obtained for this storm as It traveled 
through the Gulf of Mexico on September 22, 1975. For these data, cloud 
motion measurements were made by meteorological analysts using the Atmo- 
spheric and Oceanographic Information Processing System (AOIPS) described 
by Billingsley (1976) . These measurements have been described by Rodgers 
et al. (1979) and we have used them in evaluating the automatic computa- 
tions. As the SRI system has been discussed in considerable detail by 
Wolf et al. (1977), in this report we give only those details needed to 
describe the recent changes to it. Also, to compare our computations 
with the well-known cross-correlation method (see for example Lease et al., 
1971; Bauer, 1976; Hasler et al., 1976), we have Imbedded a cross- 
correlation computation In our methodology. This allows us to track 
cloud targets in two different ways— i.e. , using the MOTION program or 
using cross correlation. 


II CHANGES tIADE TO THE AUTOMATIC METHOD 


The major seeps in the. automatic method are shown in Figure 1. 

In the present study, substantial cimnges have been made In Steps 4 and 6. 
To separate trackable cloud groups from background (Step 4) we must first 
calculate a smoothed picture. Tills computation Is done by the TRIANG 
routine, wlilcK uses triangular weights for elements within a specified 
window (usually taken as 21 by 21 elements) around each pixel. The 
previous version, based on logic similar to that used In fast Fourier 
transforms, allowed the window to '’wrap around" at the end of a line or 
column. This distorted the smoothed picture within a distonce of one- 
half tlie window width from the edges. The new version does not allow 
wraparound, and gives a smoothed picture that is reliable to within two 
or three elements from an edge. This change aids substantially in reduc- 
ing the number of spurious motion vectors computed near edges. 

The subroutine DEVIATE in Stop 4 computen the difference between the 
brightness (or infrared value) at each point and the smoothed brightness 
at the same point. Then it selects 50 percent of the points that are 
brighter than the smoothed picture and organizes them Into "touching 
groups," A "touching group" is defined as all points that touch at least 
one other member of tlve group, with touching defined as being adjacent to 
the left, right, up, down, or along diagonals. This procedure usually 
gives some groups that are too large and irregular to serve as good 
tracers. The DELETE algorithm eliminates extraneous points from large 
groups, but for small groups (less than 10 points), all points are 
retained. For the larger groups, the proportion of points eliminated 
rises to 50 percent for groups having 50 points or more. The DELETE 
algorithm is normally repeated three times. 

We have found that the logic for eliminating points from groups is 
critically important in giving good results. In the old version of 
DELETE, a cumulative brightness distribution was made for each group. 

Then a threshold v\fas set at the particular brightness value sucll that 
the desired number of points In the group was retained. This dropped 
the dimmest members of the group, and often separated the group iiito two 
or more parts (created two smaller groups), as desired. However, it 
also had the undesirable property that it tended to discard points along 
cloud edges. The new version of DELETE operates on deviations in bright- 
ness rather than on brightness Itself. Since deviations in brightness 
tend CO be large along cloud edges (such as the edges of cirrus bands), 
the improved algorithm retains relatively bright edge targets of the sort 
tracked by human analysts. This change has added significantly to the 
accuracy of the methodology. 

To demonstrate this change, we applied DELETE to infrared data for 
hurricane Eloise. Figures 2 and 3 show visible and infrared data for an 
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FIGURE 1 BLOCK DIAGRAM OF THE AUTOMATIC CLOUD TRACKING SYSTEM 
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FIGURE 3 INFRARED IMAGE FOR THE SAME TIME AND AREA AS FIGURE 2 
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area appjroximacely 560 kro north-BOUtU by 1000 km eoBt-wust with tht» eye 
(rather IntUatinet) In the southeaBt corneri FiKuroB 4, 5, and 6 ahow 
the groupB that remain after BueeeBBive applieatlonH of DEhBTE. In 
Figure 6, the groups repreaent the cirrus along the north-Houth edge of 
the cirrufl shield and other l-fally cold cloudB, In the next step 
(Figure 7), groep» that ore very Bimill (having Iobh than two meml)crB) 
ore dlBcardcd. Also, to ovoid edge eontoinlnotion in the inotlonB, we 
eliminated groups thou touch on edge* 

If the grouping is done using vlniblo rather than infrared data, 
the results ore onalogous to tliose shown nbove. However, due to differ- 
ences between Figures 2 and 3, the visible groups are genotally not the 
same as the Infrared groups. 

As described by Wolf at al. (1977), the MOTION routine pairs groups 
on successive pictures, This is done by matching a group at time 1 with 
the best likeness of itself (in terms of else, brightness, and shape) at, 
time 2, and also requiring that the resulting motion vector conform 
approximately to the overage motion within the area being considered. 

This rationole is intended to follow the tracking concepts used by human 
analysts. The recent chonges wade In Step 6 of Figure 1 enable the MOTION 
routine to compute cloud motions simultaneously In four separate layers. 
The cloud groups are divided into four categories using the mean infrared 
value for each group based on threshold values of less than lOl, 101 to 
140, 141 to 180, and greater than 180, respectively. These categories 
tend to separate the clouds into low, middle, high, and deep convective 
types. The rationale for this very simple scheme is that the touching 
groups represent locally broken or overcast clouds} therefore the average 
infrared values are closely related to the cloud- top temperatures. 

With the groups separated into four categories, the MOTION program 
computes an average motion for each layer. •These averages normally 
differ from each other. In the iterations within MOTION, the averages 
are used to determine whether each candidate motion conforms suf ficienc.ly 
well to the associated average for that layer to he a member of the 
family. The height category is printed beside each vector in subsequent 
figures. 


Figure 8 shows the MOTIONS computed for a pair of rapid-scan pic- 
tures separated in time by approximately 10 minutes. The grouping was 
done using infrared data, and the visible data were used as an additional 
tag in the MOTION computations (see Wolf et al., 1977). The low clouds 
are traveli.ig generally southward with an average speed of 13 m a"!, and 
the high clouds (layers 3 and 4) ate moving at an average speed of 18 m s-1, 
generally toward the north. The maximum computed speed is 42 m s"l (82 
knots) for the northernmost vector in layer 3, which is for a target 
lying in the edge of the cirrus clouds. Figure 9 shows motions computed 
between a rapid-scan pair approximately one-half hour after the first 
pair. The results are quite similar. Figure 10 shows motions computed 
between two slightly later pictures separated by 30 minutes Instead of 
10 minutes. Fewer vectors are computed for the longer time interval 
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FIGURE 6 CLOUD TARGETS AFTER TWO ITERATIONS 



FIGURE 7 CLOUD TARGETS AFTER THREE ITERATIONS AND WITH GROUPS THAT 
ABUT EDGES DELETED 





FIGURE 8 CUOUD MOTION VECTORS COMPUTED FROM THE GROUPS OF FIGURE 7 

AND THE NEXT RAPID-SCAN IMAGE APPROXIMATELY 10 MINUTES LATER, 
SHOWING CIRRUS OUTFLOW AND OPPOSITE LOW-LEVEL FLOW 
Tho numbers 1 to 4 designate cloud heights. 
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FIGURE 1 1 THE SERIES OF CLOUD MOTION VECTORS COMPUTED USING SIX SUCCESSIVE 
RAPID-SCAN IMAGES EXCEPT FOR THE LAST PAIR. WHICH HAS THE NORMAL 
30-MINUTE INTERVAL 





bocauso local dynamic changos In the clouds prevent the MOTION program 
from making reliable matches for many groups. Of course the motions are 
approximataly three times larger than in Figures 8 and 9. Figure 11 
summariiics the cloud motions computed using the sec|uence of six tapes of 
8-km resolution data including the 30-minute time difference. It can be 
seen that the motions have a consistent pattern; however, there are a 
few vectors, probably erroneous, that are ni, . dds with their neighbors. 

Figure 12 shows tne imtions computed using grouping based on visible 
data for the same time as in Figure 8. Siuvc the visible data (Figure 2) 
do not have as much contrast between low and high clouds as the infrared 
data contain-, the computed field of visible motions in the center of the 
picture is rather confused and may be unreliable. For this sequence, 
better results were obtained using Infrared grouping. 

These figures illustrate the principal ImprovenKaits that hove been 
made to the SATS programs during the present contract. These improve- 
ments are crucial in obtaining accurate results in cases of complicated 
and rapidly changing clouds. For simpler iiKations, such as those associ- 
ated witli the translation of a single layer of clouds, the previous 
version of SATS performed satisfactorily. We are confident that the 
improvements are also beneficial, though to a lesser extent, for such 
cases. 



FIGURE 12 CLOUD MOTIONS COMPUTED USING VISIBLE GROUPING OF 8 km RESOLUTION 
DATA FOR THE SAME PICTURE PAIR AS IN FIGURE 8 



Ill INTRODUCTION OF A CROSS-CORRELATION COMPUTATION 


Tha most conunon method for making objective computations of cloud 
displacements uses the cross-correlation method described by Leese (1971). 
This method has been implemented for interactive use on McIDAS (Smith, 
1975) and on AOIPS (Hasler et al., 1976). To learn whether the MOTION 
program gives the same results as would be obtained by applying a cross- 
correlation computation to the same area of a picture, we have introduced 
a correlation computation as an option in SATS. The computation uses the 
groups produced by the DEVIATE and DELETE routines. These groups ore 
targets that are analogous to the cloud features selected by human 
analysts. In our usage, the target is taken as a 7 by 7 pixel area 
centered on the position of a group. The correlations are not computed 
at every point within a large area of picture 2 as is usually done. 
Instead, a search routine is used. At the start of a picture sequence, 
the search on picture 2 begins at the original location of the array on 
picture 1. Correlations of picture 1 with picture 2 are computed 
with the 7 by 7 array centered on the same point and the four neighbors 
(to the right, left, up and down). The center of the search goes to the 
point of highest correlation. Correlation computations are made again, 
and the search continues until it cannot find a higher correlation at a 
neighboring point. Then correlations are computed at all points within 
a 3 by 3 (or 5 by 5) area around the peak value. From these values, the 
point of highest correlation is found by fitting a second-degree poly- 
nomial to the nearby correlation coefficients. 

After compui-ations are made for the first pair of pictures in a 
series, the average motion in each cloud layer is computed. This motion 
is used as a first guess in the search on the next pair of images; this 
improves the accuracy of the results. The correlation computation has 
an option for using visible data, infrared data, or both. In the latter 
case, the data are simply combined since their average values and ranges 
are approximately tha same. This search type of correlation procedure 
follows that used in landmark matching by Hall et al. (1972), The com- 
putations are not made by a fast Fourier transform, but are economical 
because the search procedure restricts the number of positions that must 
be examined. In subsequent figures showing vectors determined by cross 
correlation, a result is given for each target except at edges of the 
areas. We have not used the magnitude of the correlation coefficient or 
other information in order to edit the motions. 

Figure 13 shows the motion vectors computed using the correlation 
routine for the same groups of visible data that were tracked by the 
MOTION routine in Figure 12. There are approximately twice as many 
vectors as in Figure 12 because every potential target has been used, 
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FIGURE 13 CLOUD MOTIONS COMPUTED BY CROSS CORRELATION FOR THE SAME 
GROUPS OF VISIBLE DATA USED IN FIGURE 12 


whereas the MOTION program eliminates a substantial proportion of targets 
(those that do not give acceptable motions). The pattern of motions 
appears to be reasonably satisfactory. 
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IV 


FURTHER RESUI/rS OF COMI’UTATIONS 


To invest IgaiQ the use of higher resolution data, we made cotnFuta** 
tlons for an area that comprises the left central part of Figure 3 using 
4-kin-resolutlon data and Infrared grouping. As shown in Figure 14, the 
results appear to be quite satisfactory In representing both the high- 
level and low-level clouds. For these vectors the average speeds are 
ll m s“l in layer 1, 18 m s”l in layer 3, and 15 m s-l in layer 4, There 
are 24 vectors in this area compared to 17 obtained using the coarser 
resolution data (Figure 8) . Only about 8 high-level vectors were obtained 
within the region of Figure 14 by Rodgers et al. (1979). However, if the 
resolution were increased further (e.g., from 4 km to 2 km), the auto- 
matic system would probably not be able to cope with the motions because 
there is a limitation in the initial search by MOTION — namely that the 
displacement must be less than 6 pixel units in the time interval between 
the pictures used. 


The cloud motions computed using the correlation technique for the 
same infrared groups used in Figure 14 are shown in Figure 15. The data 
coverage Is good, but the accuracy is rather marginal. In the second 
picture pair, the average layer values from the first pair are used as 
the initial guess. The rosultj! for the MOTION routine are shown in 
Figure 16, and for the COREL routlui.* in Figure 17. The latter is con- 
siderably improved compared to Figure 15. Similarly, Figures 18 and 19 
show the computations for the third picture pair. 

Figure 20 shows the cloud motion vectors computed for the series of 
six images, including a 30-minute separation of the last pair, using 4-km 
targets followed by the MOTION program. The results are generally con- 
sistent and accurate. Numerous targets are tracked over several frames. 
The comparable vectors computed by cross correlation are shown in 
Figure 21. Again the results are good, both in accuracy and coverage. 

Toward the southern periphery of Eloise the cloud patterns are 
different than to the north of the eye, being mainly cumuloform but also 
Including some high clouds. Figure 22 shows the 4-km visible image for 
this area, and Figure 23 shows the infrared data. Figure 24 shows the 
motion vectors for the first pair of pictures (separated by approximately 
10 minutes) based on automatic grouping of the visible data. There is a 
predominant strong motion of low clouds from the northwest. The average 
speed of the 18 low-cloud (layer 1) vectors is 18 m s~^-. The comparable 
motions computed by cross correlation are shown in Figure 25. The pattern 
of low cloud motions is rather confused. For the next pair of pictures, 
the layer averages of the motions from the previous picture pair were 
used as a first guess by both the MOTION and COREL routii\es. The result- 
ing motion fields are shown in Figures 26 and 27. The use of a first 
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FIGURE 14 CLOUD MOTIONS COMPUTED FOR THE LEFT CENTRAL PORTION OF FIGURE 2 
(INCLUDING THE EDGE OF THE CIRRUS SHEET! USING 4 km INFRARED DATA 
Compare to the coarser resolution of Figure 8. 
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FIGURE 15 CLOUD MOTIONS COMPUTED BY CROSS CORRELATION FOR THE SAME 
INFRARED GROUPS USED IN FIGURE 14 
















FIGURE 20 


SERIES OF SIX IMAGES OF 4 km INFRARED 
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FIGURE 21 


CWUD MOTIONS COMPUTED BY CROSS CORRELATION FOR THE 
GROUPS USED IN FIGURE 20 
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FIGURE 22 VISIBLE IMAGE USING 4-km DATA FOR AN AREA 280 by 480-km TO THE 
SOUTHWEST OF THE EYE OF ELOISE 



FIGURE 23 IMAGE OF 4 km INFRARED DATA CORRESPONDING TO FIGURE 22 
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FIGURE 27 CLOUD MOTIONS COMPUTED BY CROSS CORRELATION FOR THE SAME 
GROUPS USED IN FIGURE 26 






guess improves the results. Figure 27 has the better coverage of 
vectors. The results for the third pair of pictures are shown in 
Figures 28 and 29. Figure 30 shows the cloud motions for the five pic- 
ture pairs of the series computed by the MOTION routine using visible 
data. The patterns are quite acceptable. Figure 31 shows computations 
for the same targets computed by the COREL routine. The spurious vectors 
were computed mainly for the first picture pair. The data coverage is 
more dense than In Figure 30. Many targets can be identified in Figures 
30 and 31 as having ' een tracked quite similarly by the MOTION and COREL 
routines. The vectors computed by the MOTION routine using infrared 
targets are shown in Figure 32. There are fewer low-cloud vectors than 
in Figures 30 and 31 and more high-cloud vectors. These data, as well 
as others for cases of multilayer clouds, indicate that to obtain com- 
plete motion fields in all layers, it may be best to make tracking compu- 
tations using both infrared and visible data, if available. 

Judging the accuracy of these computed motions is difficult — as is 
almost always true for cloud motions in severe weather outside the 
radiosonde network. There are, however, approximate methods that may be 
used in this case. The first method uses a movie of Eloise prepared by 
William Shenk of NASA Goddard. By viewing this movie repeatedly or by 
using a TV display system like AOIPS or the SRI cloud console (Serebreny 
et al., 1970), one can study any sector and compare the automatically 
computed motions qualitatively to one*s impression. We have used both 
the movie and the console and feel that the results shown earlier are 
generally accurate. The second method is to compare our computations 
to cloud motions determined by Rodgers et al. (1979) for approximately 
the same period of the history of Eloise. Their computations emphasized 
visible data with l-kro and 2-km resolution and tracking by skilled 
analysts. The sequence of rapid-scan images that they used began at 
1842 GMT, approximately 45 minutes before the beginning of our sequence. 
Nevertheless, our results are in good qualitative agreement with theirs. 

For example, in Figure 8 we computed a maximum speed of a cirrus 
target of 42 m s"l (82 knots). In this same area, Rodgers et al, (1979) 
show speeds from 75 to 85 knots (see their Figure 6), In Figures 25 
through 28, the average speed of the low clouds computed with the MOTION 
routine is 11 m s”^ (21 knots) , and with the COREL routine is 10 ro s”l 
(19 knots). In this same part of Eloise, Rodgers et al. (1979, Figures 
1 and 2) show somewhat higher speeds (between 20 and 50 knots). In 
Figure 32 we show a slow drift of high clouds toward the northeast, and 
Rodgers et al. show similar values in their Figure 5. A detailed 
quantitative comparison could best be made by putting both sets of 
vectors into AOIPS and viewing them against the animated cloud sequence. 
This was beyond the scope of this study. 

In spite of these limitations in validating the automatically com- 
puted cloud motions, we believe that they appear qualitatively very 
satisfactory, and competitive in accuracy and spatial coverage to those 
obtained by Rodgers et al. (1979). 









FIGURE 32 SEQUENCE OF MOTIONS FOR FIVE PICTURE PAIRS COMf’UTED BY THE 
MOTION ROUTINE USING INFRARED DATA 


To determine to wliat extent the results of the MOTION and COREL 
routines differ, we selected a set of targets with which to compare the 
computed cloud motions. We chose targets for which both methods had 
provided satisfactory results in terms of motion vectors that were in 
approximate agreement with other nearby vectors. Also, we selected 
only those targets for which the correlation coefficient was greater 
than 0.85. The data pertain to low clouds as shown in Figure 22 
[Group (a) ] , middle and high clouds as shown in Figure 22 [Group (b) ] , 
and high clouds (including the edge of the cirrus sheet) shown in 

Figure 3 [Group (c)]. The results are given in Table 1. For Group (a), 

the average motion vectors for the MOTION and COREL routines differ by 
only 0.4 m s~l in speed and 1 degree in direction, but the rms differ- 
ence between individual vectors for the same targets is 6.7 m s“^ (0.5 
pixels in ten minutes) . For Group (b) , the average motions differ by 
2.0 m s“^ in speed and 6 degrees in directions, and the rms difference 
is 4.0 m s"l. For Group (c) , the average motions differ by 1.3 m s'"!- in 
speed and 3 degrees in direction, and the rms difference is 11.4 m s“!^. 
It would be desirable to compare these computed vectors with visual 

tracking of the same targets; however, this was beyond the scope of our 

work. We believe that the level of agreement of the averages is satis- 
factory, and that the rms differences (which arise from differences in 
location of only fractions of a pixel) represent largely an unavoidable 
random component in computations of this sort. 
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Table 1 


COMPARISON OP CLOUD MOTIONS (units m s"!) 
COMPUTED USING THE MOTION AND CORRELATION ROUTINES 


(a) Low Clouds, 25 Pairs, 4-km Visible Data (Figure 22) 

Average 

Motion 

Correlation 

RMS Difference 

u 

15.5 

15.8 

3.9 

V 

-19.6 

-18.8 

5.4 

Speed 

24.9 

24.5 

6.7 (0.5 pixels) 

Direction 

321“ 

320“ 


(b) High Clouds, 19 Pairs, 4-km Visible Data (Figure 22) 

Average 

Motion 

Correlation 

RMS Difference 

u 

9.1 

11.4 

3.6 

V 

-3.8 

-3.2 

1.9 

Speed 

9.8 

11.8 

4.0 (0.3 pixels) 

Direction 

292* 

286“ 


(c) High Clouds, 42 Pairs, 4-km III Data (from Figure 3) 

Average 

Motion 

Correlation 

RMS Difference 

u 

-7.1 

-8.4 

8.8 

V 

30.4 

28.7 

7.1 

Speed 

31.2 

29.9 

11.4 (0.8 pixels) 

Direction 

167“ 

164“ 
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V SPMMAKY ANI> i t^Nri.US IONS 


111 thin roHiMrrh, tin* SKI .itil i>m.i I it* rliMul ti.ii'kiiif* im^tluulM wm* 
siil>H I ant l.i I I V i mptiWiHl • Altlioii>»li t lu* trsls m.uli> an* llmiti*il in ntimlu*r« 
tin* 4*ast»H i'htts(*n 1 4»r liiir r I i*aiu» Kli>lH4* .in* V4*rv 4litlU*nlt oii 4 *m 4*V4*n 1 4*t 
human ir.u'kinK; iu*V4»r t Ii 4 * I 4»hh , I lu* anlomati4* HVHn*m |u*r 1 4*rm4*4l V4*rv W4*ll. 

In simpU* r.iS4*s tli4*t4* *ipp4*.ii*H to lu* llttU* 4liMiht 4>t tin* a4*4*tira4*y 4M tin* 
4 * 4 Mtiptit 4*4! 4 ‘liuiil moti 4 MiH, .iM W 4 * r 4 »port 4 * 4 l 1 4 ii tin* pr4*viiuiH V4 *ikI 4MI 4>l tin* 
m4*t lnul4> lo>*v (W4»lt 4»t . 1 1 • • hl/y). In 4umip I i 4*.it 4*4! 4MS4*h, H4mu* 4UMitr4>l 
oV4*i or 4*4litinp, o\ aiit imiat i4*.i I ly lunnpiit 4*4! mi>t hnis mav l> 4 * iu* 4 * 4 l 4 * 4 l • 

rin* lu*st r4*siiltH in t4*rms 4>! .uu'iii .n*v aiul 4*4>V4»ra>»4» tor tin* liuri ii'.iin* 
riols4* 4lata W4*i 4 * 4>ht.iiin*4l iisin>; /»-km-r4*Si»liit I 4 MI ilat.i at I O-ml niit 4 * Int4*r- 
vals. Ki»i pli*tiii4*s S 4 *p.ir.it 4 * 4 ! I>v lnt 4 *rv.ilM 4)1 approximat 4 *ly 10 ml»nit 4 *s, 
tin* aiitomati4* 4*4)mpnt at i 4)iis r4*4pili4* 4lat«i with 4*o.irs4*r r4*solnt l4>n in 4)ril4*r 
to 1)4* ahli* t4) t r.i4 k 1 * 4 * I .it i V4* I v h i ^Ii-sp4*4*4l nu>t ions. In 4)tln*r wor4ls« tin* 
nu*thoii han tia* limit.it i4)ii tli.it t inu* .iiul sp.nu* S4*al4*s must l>4* m.it4*ln*4l S4) 
tli.it tin* I'loiiil 4 I i sp I.i 4 * 4 *m 4 *nt s .in* U*sh th.in t> pix 4 *l units l) 4 »tW 4 * 4 *n pi 4 *tiir 4 *s; 
it in)t , .in appri'x imat 4 * tirst ^iu*ss 4)t tin* il ispla4*4*m4*nt is iu*4*il4*il si) th.it 
tin* Mi>rii>N ri)utiiu» or tin* (HtKT.h routiiu* will S4*.ir4*h in tin* ripjit l4)4*.ili- 
l i4»s 1 4 ) 1 ’ m.it4*hinp« >;r4»iips. I\) .1 i*4*rt.iin 4»xt4»nt in sinMi i*.iS4*s, SAI’S 
impri)V4*s its tirst >*U4*ss t r4»m tin* tirst 4)in* 4 ) 1 * tW4> pi4*tur4* p.iirs 4>t «i 
S4*4jU4*iu*4» , S4) tli.it tr.u'kinp, is r4*li.ihl4* in l.iti*r p.iirs. A 1 1 4»rn.it i V4* 1 v , 

4'4>r put .It i4)iis m.ul4* with i*4>.irs4*-r4*S4) hit i 4)11 4l.it.i i*.in l>4* us4*il .is .1 tirst 
^ii4*ss 1 4»r nu'r4* il4*t .1 i h*4l .iiui .i4*4*urat4* 4*4»mput .it i4»ns with tiiu*r r4*S4)lut i4>n 
4 l.it .1 • 

An impi'rt.int new opt K^n h.is hoen .ukleJ ti' the methi)ili>l 4 >gv in terms 
4)f the cross-4.'4>rre lat ion computation i)f moti 4 >ns. In this 4 >pti 4 >n the 
t.irgets are selecteil autimiat ical ly, but the displacements are ci^mputed 
hv .1 se.irch r4>utine that U)cates the pi'int of highest C 4 >rre I at ion- - i , e , , 
the OOKKL nniline replaces the MlvriON routine. Both rout ines give satis- 
factory results. For typical targets the displacements computed by the 
NOTION and CORKL routines differ by luily fractions of a pixel. In our 
opinion this level of accuracy is as high as one can expect since the 
data contain sm.ill-scale time-varying features. 

At this time we believe that more extensive tests should be made by 
including SATS in an interactive system such as AOlPS, where it can be 
applied to a large number of cases. The results can then be checked 
directly against measurements made under operati'ir ci^ntrol. Experienced 
analysts should then be able to make a comprehensive Ci>mparison and 
evaluation of the different options and methods. 
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